Quasiparticle-scattering measurements of laminar and turbulent vortex flow in the 

spin-down of superfluid 3 He-B 
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The dynamics of quantized vortices is studied in superfluid 3 He-B after a rapid stop of rotation. 
We use Andreev reflection of thermal excitations to monitor vortex motion with quartz tuning fork 
oscillators in two different experimental setups at temperatures below 0.2T C . Deviations from ideal 
cylindrical symmetry in the flow environment cause the early decay to become turbulent. This is 
identified from a rapid initial overshoot in the vortex density above the value before the spin-down 
and its subsequent decay with a t~ 3 ^ 2 time dependence. The high polarization of the vortices along 
the rotation axis significantly suppresses the effective turbulent kinematic viscosity v' below the 
values reported for more homogeneous turbulence and leads to a laminar late-time response. The 
laminar dissipation down to T < 0.15T C is determined from the precession frequency of the polarized 
vortex configuration. In the limit of vanishing normal component density, it is found to approach a 
temperature-independent value, whose origin is currently under discussion. 

PACS numbers: 67.30. hb, 47.15.ki, 47.20.-k, 67.25.dk 
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I. INTRODUCTION 

The dynamics of quantized vortices, especially super- 
fluid turbulence at ultralow temperatures, where the vor- 
tices can move with little dissipation, is an active topic of 
research^ In the T — > limit, the superfluid can be mod- 
eled as an inviscid and incompressible fluid, where turbu- 
lence involves vortex reconnections and tangle formation. 
In these conditions, numerical modeling and the inter- 
pretation of measurements on turbulence become more 
manageable than in the case of classical viscous fluids. 

One way to generate complex fluid motion is to bring 
a rotating container to rest. In classical fluids, the re- 
sponse of the fluid in such a spin-down experiment is 
governed by the frictional forces at the boundaries and 
becomes unstable at relatively low Reynolds numbers^ 
Also in superfluid 4 He, where vortex pinning at the con- 
tainer walls plays a significant role due to the small vortex 
core radius, an impulsive spin-down to rest is generally 
expected to generate turbulence^ 

In superfluid 3 He the situation is different, since in the 
case of smooth walls, the coupling to the container is ac- 
complished mainly by the frictional interaction between 
the normal and superfluid component known as mutual 
friction. That is, the superfluid does not couple to the 
boundaries directly, but instead, does so through vol- 
ume forces to the normal fluid, which is stationary in the 
reference frame of the container. At high temperatures 
and large mutual friction vortex motion is laminar, but 
as the mutual-friction coupling vanishes exponentially in 
the zero-temperature limit, the flow of vortices is easily 
destabilized and becomes turbulent. Earlier studies in 
3 Hc-B show that in an axially symmetric environment 
the spin-down flow is laminar even in the limit of vanish- 
ing normal fluid density— However, with increasing sur- 
face friction, the response ultimately becomes turbulent^ 
Here we show that a similar change from a laminar to 



a turbulent spin-down response can be brought about 
by deviations in the flow geometry from axial symme- 
try. The ensuing turbulence is concentrated in the initial 
phase of the response and since the polarization of the 
vortices along the symmetry axis remains high, laminar 
flow ultimately takes over the late part of the response at 
lower vortex densities even at temperatures below 0.15T C . 



In superfluid 3 He-B, vibrating objects can be used to 
monitor the density of thermal excitations. This has 
proven to be useful for studying different vortex struc- 
tures, especially at ultralow temperatures, where prac- 
tically all other detection methods become insensitive.— 
The measurements rely on the fact that the rotating flow 
associated with quantized vortices can constrain the tra- 
jectories of the elementary excitations, which otherwise 
move along ballistic flight paths. In the presence of su- 
perfluid flow with velocity v s excitations with momentum 
p undergo the Galilean transformation E — > E + p • v s . 
An excitation moving in the superfluid with insufficient 
energy to overcome the potential energy barrier created 
by the flow of vortices has no forward-propagating states 
and therefore retraces its trajectory changing flavor from 
quasiparticle to quasihole in a process called Andreev 
reflection^ Thus, the presence of vortices can be inferred 
from the variations in thermal damping of vibrating ob- 
jects, such as vibrating wires or quartz tuning forks. The 
cross section of Andreev reflection is large even for sparse 
vortex structures providing a sensitive probe of the super- 
fluid flow fieldi£i£ Here we use these techniques to study 
the interplay of laminar and turbulent flow of vortices in 
the spin-down response when a deviation is introduced 
from ideal axially symmetric flow. 
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II. TURBULENT AND LAMINAR DYNAMICS 
OF SUPERFLOW 

A. Turbulent superflow 

Turbulence is often defined as a complex and dynamic 
flow field, which involves processes spanning several or- 
ders of magnitude in spatial extent with aperiodic tempo- 
ral dependence. In quantum turbulence, different types 
of turbulent flows are characterized by the local density 
L(r) of the vortices and their polarization. In its simplest 
form, turbulence in superfluids consists of a homogeneous 
and isotropic tangle of singly quantized vortex lines. 

On length scales large compared with the intervortex 
distance £ = L -1 / 2 quantum turbulence often resembles 
its classical counterpart* 1 - The energy is injected into ed- 
dies at length scales determined by the characteristic 
size of the flow disturbance. The large scale motion is 
achieved by partial polarization of vortices to bundles 
forming eddies of different sizes. As in classical turbu- 
lence, the energy cascades down with a Kolmogorov-type 
energy spectrum given by 



E(k) 



Ce 2 / 3 fc 5 / 3 , 



(1) 



where C sa 1.5 is the Kolmogorov constant. Assuming 
the dissipation is determined by the length scale £, the 
energy flux towards shorter length scales, i.e., inverse k, 
per unit mass is given by 



/ 2 T i 
e = v k L 



(2) 



where v' is the effective kinematic viscosity and k the 
circulation quantum. Using Eqs. ([J) and ([2]) allows us 
to describe the late-time decay of vortex density^ by 



27TV V 1 'K 



-3/2 



(3) 



where the container size D determines the cutoff 
wavenumber fco = 2-k / D. The value of v' depends both 
on temperature and the nature of the flow. Despite the 
fact that the derivation leading to Eq. © assumes homo- 
geneous and isotropic turbulence, our experiments show 
that turbulence with high polarization of vortices can 
also decay as i~ 3 / 2 . However, the results suggest that 
polarization of the turbulent vortex structure suppresses 
v' significantly compared to the value measured for more 
isotropic turbulence. 

It is widely believed that as the energy flows to length 
scales smaller than the intervortex distance, it is trans- 
ferred in a cascade of helical deformations of individ- 
ual vortices called Kelvin waves X There is an ongo- 
ing debate on the nature of the energy transfer from 
the classical Kolmogorov-like cascade to the Kelvin- wave 
cascadefii - — At finite temperatures mutual friction pro- 
vides dissipation at all length scales, and even at T = 0, 
the energy is ultimately dissipated when Kelvin waves at 
very large k induce phonon emission^ in 4 He or perhaps 
quasiparticlc emission^ from vortex cores in 3 Hc-B. 



B. Laminar superflow 

The other extreme of dynamic vortex motion is fully 
laminar polarized flow supporting no vortex reconnec- 
tions. In supcrfluid 3 Hc-B the normal component is 
clamped to corotation with the cylindrical container with 
velocity v n = ilzxr in the cylindrical coordinates (r, 0, z) 
with 51 being the angular rotation velocity. At ultralow 
temperatures, where the reactive mutual friction can be 
neglected, the course-grained hydrodynamical equation 
for the superfluid velocity v s is given byi£ 



dt 



+V^-(v s -V)v s = -awx((v s -v n )x(Vxv s )), (4) 



where fi is the chemical potential and o> a unit vector 
along the vorticity V x v s . In the ballistic regime, the 
dissipative mutual-friction parameter a is proportional 
to the Boltzmann factor exp(— A/T), where A is the 
pressure-dependent supcrfluid energy gap. 

If the vortices remain highly polarized along the axis of 
the cylinder, and thus vortex reconnections play no role, 
laminar rotating flow of the supcrfluid is solid-body-likc 
with velocity v s = O s z x r. In this case, taking a curl of 
both sides simplifies Eq. © to 



dft s (t) 
dt 



2aCt(t)p(t) - n s (t)} 



(5) 



For a step change of the rotation drive at t = from 
an angular velocity fio to rest the solution of Eq. <j5j> is 
given by 



n.(t) 



n 



i + t/i 



(0) 



where r = (2ai7o) _1 . In real experiments the change 
of rotation velocity is done at a finite rate, in our case 
typically at a = dVt/dt ~ —0.03 rad/s 2 . During the de- 
celeration, i.e., for — ilo/a < t < 0, the solution of Eq. 
© is 



n.(t) = 



,a(t+fi /a)(fi -at) 



To\/d + y / 7rae Qfi o/ a crf( v / aat) 



(7) 



where tq = Qq 1 + y/a/aexp(ail,Q/a)eTi(y/a/aQo). At 
low temperatures, where a <C a/fip the supcrfluid veloc- 
ity at the end of the deceleration is very close to £Iq and 
consequently, Eq. © can be used at all times t > 0. 



III. MEASUREMENT TECHNIQUES 

In our present measurements the 3 He-B sample is con- 
tained in a smooth-walled fused quartz cylinder. Two 
experimental setups shown in Fig. [T]have been used. In 
setup 1 on the left the 3 He-B sample is pressurized to 
29 bar and the cylinder is divided in three parts by two 
0.7-mm-thick quartz division plates with 0.75 mm and 
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FIG. 1. (Color online) (Left) Experimental setup for the 
bolometric measurement of the early turbulent spin-down 
response. The upper experimental volume modeled as a 
bolometer is separated from the heat exchanger volume at 
the bottom by the lower division plate with a small 0.3-mm- 
diameter orifice. The upper division plate with a 0.75-mm- 
diameter aperture secures laminar spin-down flow in the top- 
most section. The quartz tuning forks are calibrated to mea- 
sure the heat transport to the refrigerator. (Right) Setup 
without division plates for the measurement of the late lami- 
nar part of the spin-down response. 



0.3 mm orifices. The lowermost part consists of a 30- 
mm-long, 3.6-mm-inner-diameter tube, which opens at 
the bottom to the heat exchanger. This volume is cooled 
to the lowest possible temperature T < 0.14T C with a 
nuclear demagnetization cooling stage. The two upper 
parts combined constitute a 12-cm-long smooth-walled 
section of the quartz tube with 6 mm inner diameter and 
act here as a bolometer. The top section above the orifice 
of 0.75 mm diameter has been used in earlier measure- 
ments for examining the stability of laminar vortex flow 
with NMR techniques^ The middle section located be- 
tween the two division plates was added later to measure 
Andrccv reflection from a cluster of rectilinear vorticcsA£ 
and the thermal signal from vortex dissipation in turbu- 
lent spin-up.— Here we use the tuning-fork oscillators in 
the middle section for measuring the spin-down response 
in the bottommost section. In this case vortex flow is 
perturbed at the bottom end of the sample tube by the 
presence of a rough and wavy sintered heat-exchanger 
surface and the leads of the tuning forks. 

In setup 2, shown on the right in Fig. [1] both divi- 
sion plates are removed and the sample is pressurized to 
0.5 bar instead of 29 bar. In this setup, axial flow symme- 
try is disturbed by the presence of the two tuning forks 
placed at the bottom end of the quartz tube. 



In both setups the container is initially filled with a 
uniform array of rectilinear vortices oriented along the 
rotation axis. To create the vortex array the cryostat is 
rotated at constant angular velocity H around the axis 
of the container tubeJ^ The vortex density in the equi- 
librium state is determined by minimization of the free 
energy in the rotating frame and is given by the solid- 
body-rotation value L — 2Q/k. 

In the absence of applied heating, the temperature 
T sa 0.20T C in the bolometer of setup 1 is determined 
by the background heat leak from the container walls 
(Qhi ~ 10 — 20 pW depending on rotation velocity) and 
the thermal resistance Rt(T) = (dQ / dT)^ 1 from the 
tiny 0.3-mm orifice in the lower division plate. The ther- 
mal balance in the system is discussed in more detail in 
Ref. 17]. Assuming thermal equilibrium in the volume 



above the 0.3-mm orifice, all of the heat flux Q must leave 
through the orifice as a flux of energy-carrying excitations 
with the temperature dependence given by 



Q(T) 



ink B p 2 F 
h 3 



Te-r(A + k B T)A h . 



(8) 



Here is the so-called effective area of the orifice, which 
depends on the number and orientation of the vortices 
below the 0.3 mm orifice. The reason for this is that due 
to the vortices, part of the thermal excitations experience 
Andrccv reflection and thus, the effective area is reduced 
(see the inset of Fig. [2]). Therefore, information on vortex 
dynamics can be obtained from the temporal dependence 
of A h . 

The principle of the measurement is to study how the 
Andreev scattering from the vortices affects the fraction 
of thermal excitations returning to the bolometer, 



fr(t) = l-A h (t)/A hfi , 



(9) 



where A^.o is the effective area in the absence of vor- 
tices. The temperature is obtained from the resonance 
lincwidth A/ of a quartz tuning fork oscillator, which is 
proportional to the Boltzmann factor— exp(— A/T) af- 
ter subtracting the intrinsic linewidth A/; nt « 15 mHz. 
Thus, for known Q(T), Ah(t) can be calculated directly 
from Eq. ||5J). The background heat leak Qhi in the 
steady state can be determined by measuring the tem- 
perature increase with one tuning fork as a function of a 
known heat input from the other forkjii 



IV. BOLOMETRIC MEASUREMENT OF 
TURBULENCE IN THE SPIN-DOWN RESPONSE 

Figure [5] shows the temperature response of the 
bolometer to a rapid change of rotation. The moment 
t = corresponds to the time when the final constant 
value = is attained after a change of drive at fl = 
—0.03 rad/s 2 from the initial velocity f2i n ; = 0.5 rad/s. 
The initial temperature increase during the deceleration 
is a consequence of three different phenomena: direct 
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FIG. 2. (Color online) Temperature T in the bolometer and 
an estimate for the lower limit of the vortex density L after 
bringing the container to rest from f2i n i=0.5 rad/s at t = 0. 
The vortex density is inferred from the fraction of Andreev- 
reflected thermal excitations as discussed in the text. The 
solid line is a fit to L oc (t + r) _3//2 dependence after the 
initial overshoot in the data. The inset shows the principle 
of the measurement: The vortices below the orifice reflect 
some fraction of the thermal excitations back to the bolome- 
ter. The fraction depends on the density and configuration 
of these vortices at the lower temperature T < 0.14T C . The 
measurement is performed at 29 bar liquid 3 He pressure. 



heating from the decay of vortices inside the bolometer 
volume, weakly increased rotation-induced heat leak at 
certain velocities of mechanical resonances of the cryo- 
stat in the range < fl < f2ini, and most importantly, 
increased thermal resistance Rt due to the spin-down- 
induccd turbulence below the bolometer volume. 

After the container is brought to rest, the heat gener- 
ated in the bolometer Q gen (T), which is at T sa 0.20T C , 
is the sum of the known background heat leak and the 
heat release from the laminar decay of vortices stabilized 
by the mutual friction a ps 0.002 in the upper sections 
of the cylinder above the 0.3 mm orifice.— The kinetic 
energy of the rotating superfluid with density p s is given 
by 



£ki„ = np s R A hQ 2 Ji, 



(10) 



where R and h are the radius and the height of the cylin- 
drical bolometer volume, respectively. Combining Eqs. 
© and (fTTH) yields the rate at which the kinetic energy 
translates into heat, 



^ kin = ![^M (1 + , /r) -3. 

At 



(11) 



Now the total heat flow out of the bolometer volume is 
given by 



where V = wR 2 h is the volume of the bolometer and the 
heat capacity c(T) is given approximately bySI 



c(T) 



k 2 



2ttN p 



2 1' 
16" 



-T 



(13) 



Here, Np is the density of states at the Fermi level. From 
the measured temperature and its time derivative T, one 
can solve for the reflection coefficient f r (t) using Eqs. (jHJ, 
(El), (32), and (Q3). 

Even though we do not know the detailed structure 
of the vortex configuration below the 0.3 mm orifice re- 
flecting a part of the thermal excitation beam back to the 
bolometer, two limiting cases can be considered to obtain 
upper and lower limits for the vortex density L. Obvi- 
ously, the fraction of the Andreev-reflected beam at fixed 
L is minimized if all of the vortices are polarized along the 
axis of rotation and are thus perpendicular to the plane 
of the bolometer orifice. Accordingly, the upper limit for 
the vortex density at a given reflection coefficient can be 
obtained from the measurements of / r (f2) with known 
vortex density in uniform rotation (as measured in Rcf. 
[l7j]). The lower limit can be estimated by considering a 
system of vortices which are perpendicular to the beam 
of excitations and have an intervortex spacing £ = L~ 1//2 . 
In this case the vortex density is given byS 



L = - 



2k B T\n(l-f t ) 
37rA£S* 



(14) 



Q(t) = Qgen(i) ~ C(T)VT, 



(12) 



where £ — 2 ■ 10 -8 m is the coherence length and S the 
spatial extent of vorticity. This simple estimate does not 
account for the diffusively scattering boundaries of the 
container. Our numerical simulations^ show that in the 
presence of vortices the wall effects increase / r (f2) by al- 
most exactly a factor of two. Thus, a reasonable lower 
limit for L is attained by dividing the measured reflection 
coefficient by two and taking S to be the mean distance 
for excitations to travel before they would hit the bound- 
ary of the container in the absence of vortices. Figured] 
shows an example of the lower limit of L calculated this 
way from the measured temperature response shown in 
the same figure. The upper limit obtained from the cal- 
ibration measurements of straight vortices with known 
density is a few times larger. 

The result for L(t) in Fig. [^carries the following two 
characteristic signatures from turbulence. First, an ini- 
tial overshoot, of the order of a few tens of percents with 
the maximum at around t = 30 s, indicates that part of 
the kinetic energy of the superfluid is converted to a tur- 
bulent tangle of vortices. The upper and lower limits of 
the maximum vortex density as a function of initial rota- 
tion velocity are shown in the right panel of Fig. |3] The 
lower limit maximum lies close to the known initial den- 
sity 2f2 ini /ft. Second, after a couple hundred seconds the 
estimated vortex density (both upper and lower limit) fits 
well with a t -3 / 2 dependence on time. This suggests that 
at least some fraction of the vortices decays in a turbu- 
lent manner, since for a fully laminar response the decay 
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FIG. 3. (Color online) (Left) Lower limit of the vortex density L as a function of time after bringing the container to rest for 
four initial velocities fii n i (see text for details). A simple moving average filtering is used to reduce the noise in the data. The 
upper limit of L at each value of Q.^ obtained from a calibration measurement with rectilinear vortices at a known density is 
a few times higher. (Right) Lower and upper limits for the maximum vortex density generated in spin-down together with the 
initial steady-state value 2fii n i/ft. 



should be proportional to i -1 , as given by Eq. ([5]). In 
the left panel of Fig. [31 examples of the lower limit of the 
vortex density L at different initial rotation velocities are 
plotted with logarithmic scales illustrating the i -3 / 2 time 
dependence. After ~ 10 3 s the vortex densities become 
too low to be detected with good accuracy. 

The overshoot in Fig. [5] is roughly an order of mag- 
nitude smaller than in the more turbulent spin-down re- 
sponse of the measurements in supcrfluid 4 He in a cubic 
container with grid-covered surfaces in Ref. [H, but of 
similar magnitude as in the numerically calculated ex- 
amples of more polarized turbulence in Refs. [|| and 
[j| . The characteristic time to reach the maximum of the 
overshoot in L(t) decreases with increasing initial rota- 
tion velocity £7; n ;. With the present deceleration rate, 
the time to bring the cryostat to rest exceeds the time 
to reach the maximum level of L(t) at high velocities^ 
In this situation, the line density remains approximately 
constant during the remaining deceleration. 

It is tempting to relate the response in Fig. [3] to 
Eq. §5§ for the quasiclassical decay of the line density 
with the same i~ 3//2 time dependence. One needs to 
be careful though, since some of the assumptions lead- 
ing to Eq. (J3J, especially the homogeneity of the flow, 
arc not met. Nonetheless, if the size of the energy con- 
taining length scale is chosen to be the diameter of the 
container, i.e. D = 2R, then the effective kinematic vis- 
cosity v' extracted from the lower limit of vortex densities 
is (1.0 ± 0.3) • 10~ 4 k and does not depend on the initial 
rotation velocity in the range 0.2 < fi^i < 1.7 rad/s. 
The estimate obtained from the upper limit is roughly 
an order of magnitude smaller. 

Earlier experiments on the free decay of turbulence in 
3 Hc-B in the T = limit have been done in a series of 
measurements by Bradley et al. They used a vibrating 



grid to create a tangle of vortices and inferred the vor- 
tex density from a reduction in the damping of vibrating 
wire resonators caused by Andreev reflection from the 
flow field of the vortices.— In their experiment, the decay 
was more than an order of magnitude faster and yielded 
v' w 0.3k, which is two orders of magnitude larger than 
the value measured in 4 He for the free decay of tangles 
produced by rapid spin-down.— The orders-of-magnitudc- 
smaller value for the effective kinematic viscosity in our 
experiment underlines the influence of the large polariza- 
tion of vortices, which suppresses the amount of vortex 
reconnections and thus dissipation in comparison to more 
homogeneous turbulence. 

Since the reflection coefficient is more sensitive to the 
vortices which arc perpendicular to the excitation beam 
than to those parallel to it, more quantitative estimates 
of L cannot be done without detailed information on the 
polarization of the vortex structure. Supposedly, in the 
beginning the vortex configuration is more randomized 
and less polarized along the rotation axis. These vortices 
reconnect with larger probability than the polarized vor- 
tices and account for the observed £ -3 / 2 dependence in 
the early part of the decay. This assumption is supported 
both by numerical 2 ^ and experimental (Sec. |Vj results. 

V. LAMINAR DECAY INFERRED FROM A 
PRECESSING VORTEX CLUSTER 

In the absence of the division plates in setup 2, the 
quartz tuning-fork resonators can be used to probe the 
spin-down-induced flow field locally. In the equilibrium 
vortex state at constant rotation, the ballistic nature of 
thermal excitations assures thermal equilibrium inside 
the container. Dynamic inhomogencous vortex struc- 
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FIG. 4. (Color online) Response of the resonance width 
of the quartz tuning fork oscillator to a rapid spin-down from 
f2ini=1.02 rad/s to rest. The initial increase arises mainly ow- 
ing to Andreev reflection from the rapidly decaying turbulent 
vortex configuration generated in the early part of the spin- 
down. The inset shows the zoomed view of the oscillations 
in the laminar late response caused by a periodic variation 
of the thermal excitation density in the vicinity of the fork. 
These oscillations result from the Andreev scattering from 
an asymmetric vortex cluster precessing at angular velocity 
Q. B (t). The measurement is performed at 0.5 bar liquid 3 He 
pressure. 
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FIG. 5. (Color online) Superfluid angular velocity Q s as a 
function of time after bringing the container to rest at t = 0. 
The response fits perfectly to the laminar flow model [Eq. 
©] shown by the dashed line with Qo ~ 0.80Sli n i and a ~ 
6.45 • 10~ 4 . The inset shows the dissipative mutual-friction 
parameter a as a function of the linewidth of the tuning-fork 
oscillator (bottom axis) or temperature (top axis) . The value 
of a, whose uncertainty is of the order of ±5 ■ 10 -5 , is ex- 
tracted from the precession frequency of the vortex cluster, as 
discussed in the text, and averaged over 2 — 4 measurements at 
different rotation velocities in the range fimi = 0.6 — 1.5 rad/s 
at the given temperature. 



tures, however, can shadow a part of the background 
excitation flow emanating from the walls, and thus, pro- 
duce local temperature variations. Recently, these varia- 
tions, measured with vibrating wires, were used to study 
line-density fluctuations in a turbulent vortex tangled 

Here we study the response of the superfluid, origi- 
nally rotating at constant angular velocity Q s = Oi n i, 
to a rapid spin-down to rest by following the variation in 
the thermal excitation density in the vicinity of a tuning- 
fork oscillator. The background temperature in the range 
150 — 190 /j,K is controlled by magnetization and demag- 
netization of the nuclear cooling stage. 

Figure 2] shows an example of the fork-oscillator re- 
sponse before, during, and after bringing the container to 
rest. The initial increase of the linewidth arises mainly 
from Andreev scattering from the turbulent vortex con- 
figuration generated in spin-down as described in Sec. 
HVl After about 200 s, periodic oscillations can be ob- 
served on top of the relaxing average temperature. We 
interpret these oscillations to originate from a cluster of 
predominantly straight vortices precessing at the super- 
fluid angular velocity Q s . Other experimental signals of 
a similar type from preces sing vortex bundles have been 
reported in Rcfs. |18f and [261 ]. Here a small asymmetry 
in the structure of the cluster causes periodic variation 
of the thermal excitation density in the vicinity of the 



tuning fork oscillator. The frequency of these oscilla- 
tions drops as < _1 for several hours indicating that the 
cluster decays in a laminar manner. Using rio and r as 
fitting parameters produces an excellent fit to Eq. © 
with £7o = bQini and r oc fl^ 1 as expected for solid-body- 
like laminar decay, illustrated in the right panel of Fig. 
|U Measurements at different initial rotation velocities 
^ini = 0.6— 1.5 rad/s and temperatures T = 0.15— 0.19 T c 
show that the constant b 0.8 and depends neither on 
Q- m i nor temperature. 

Both fitting parameters provide valuable information 
on the vortex dynamics during the spin-down. The 
rapid drop of the superfluid angular velocity from f2j n j to 
~ 0.8fi; n i indicates that during the turbulent early part 
of the decay about one fifth of the vortices are annihilated 
due to vortex reconnections. The initial turbulent part 
seems to last a shorter time than in the measurements 
at 29 bar in Fig. [3] Two reasons for the increased tur- 
bulence in the high-pressure measurements can be given: 
At 0.5 bar pressure the vortex core radius is roughly five 
times larger than at 29 bar. Therefore, vortex pinning on 
the surface of the heat exchanger, which is sintered from 
~ 10-/xm-size copper flakes, is expected to be stronger at 
high pressures. Second, the short narrow-diameter bot- 
tom part of the cylinder in setup 1 is markedly different 
from that of setup 2. 

The decay of the remaining polarized vortices is dom- 
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inated by the mutual friction with the parameter a = 
(2tQq)~ 1 . To date, the parameter a at low pressures 
has been measured only at temperatures above 0.35 T c ,— 
where the normal fluid density is three orders of magni- 
tude higher than at the lowest temperature of this work. 
Eltsov et al.— used NMR-techniques to measure a down 
to 0.20 T c at high pressures. 

The vortex mutual friction is a consequence of the 
interaction between the thermal and vortex-core-bound 
excitations^ Since the normal-fluid density in the bal- 
listic regime drops exponentially, the mutual-friction pa- 
rameter should be proportional to the linewidth of the 
tuning-fork oscillator, A/ oc exp(— A/T). The inset of 
Fig. [5] demonstrates how a follows the expected ex- 
ponential dependence on temperature, but with a non- 
zero intercept a(0) ~ 5 • 10 -4 in the T — >• limit. 
The source of this zero-temperature dissipation is un- 
clear. One possibility is surface interactions with the con- 
tainer walls. Another recently suggested mechanism is lo- 
cal heating of the vortex cores in accelerating motion^ 
which can lead to temperature-independent dissipation in 
the zero-temperature limit. Finite zero-temperature dis- 
sipation has been observed earlier in the turbulent front 
propagation^ Whether the origin of the finite value for 
a is the same here as in the front motion is uncertain and 
remains to be clarified. 



VI. CONCLUSIONS 

In conclusion, we have studied the spin-down response 
of supcrfluid 3 Hc-B in a cylindrical container at very low 
temperatures, where the normal fluid excitations are ex- 
tremely dilute. The low mutual friction together with 
small deviations from an ideal axially symmetric flow en- 
vironment results in turbulence. This is in contrast to our 
earlier measurements at higher temperatures T > 0.2T C 
with more symmetric geometry^ 



Vortex reconnections concentrate in the early part of 
the decay where they cause loss of polarization, increased 
dissipation, and i~ 3 / 2 time dependence of the vortex den- 
sity. Similarly to increased surface friction in Ref. 
the increased vortex pinning on the rough and wavy heat- 
exchanger surface leads to increased turbulence also. The 
vortex core size seems to affect the decay, which agrees 
with the conclusion that vortex pinning plays a role in 
maintaining turbulent flow. 

The vortex density in the late decay can be inferred 
with high accuracy from the local temperature variations 
caused by Andreev reflection from the processing slightly 
asymmetric cluster of quantized vortex lines. This novel 
method reveals that the late decay of the vortex cluster 
is laminar even at the lowest temperatures. The tem- 
poral dependence of the spin-down decay allows us to 
determine both the amount of vortices lost in the ini- 
tial turbulent part of the decay and the vortex mutual 
friction down to much lower normal fluid densities than 
before. 

The dissipative mutual friction parameter a is found 
to approach a small nonzero value in the low temperature 
limit. This is the first time this type of residual friction 
is observed for laminar motion. Whether it signals true 
zero-temperature dissipation or some nonexponentially 
decaying dissipative phenomenon remains as an interest- 
ing challenge for the future. 
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